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Abstract 

In this paper we investigate numerically the model for pedestrian traffic proposed in [B. Andreianov, C. Donadello, 
M.D. Rosini, Crowd dynamics and conservation laws with nonlocal constraints and capacity drop. Mathematical Models 
and Methods in Applied Sciences 24 (13) (2014) 2685-2722]. We prove the convergence of a scheme based on a constraint 
finite volume method and validate it with an explicit solution obtained in the above reference. We then perform ad hoc 
simulations to qualitatively validate the model under consideration by proving its ability to reproduce typical phenomena 
at the bottlenecks, such as Faster Is Slower effect and the Braess’ paradox. 
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1. Introduction 

Andreianov, Donadello and Rosini developed in |[ll] a macroscopic model, called here ADR, aiming at describing the 
behaviour of pedestrians at bottlenecks. The model is given by the Cauchy problem for a scalar hyperbolic conservation 
law in one space dimension with non-local point constraint of the form 


d,p + djip) = 0 

(C ;ii:) 6 M+ X M, 

(la) 

p(0, x) = p(x) 

X c M., 

(lb) 

f(p(t,0±))<p^J w(x) p(f, x) dxj 

t £ 

(Ic) 


where p(t, x) e [0, R] is the (mean) density of pedestrians in e M at time f e R+ and p: M ^ [0, R] is the initial (mean) 
density, with R > 0 being the maximal density. Then, /: [0, R] —> M+ is the flow considered to be bell-shaped, which is 
an assumption commonly used in crowd dynamics. A typical example of such flow is the so-called Lighthill-Whitham- 
Richards (LWR) flux ||l[l|l defined by 

P 


f(p) =pVtnax 1 


where Vmax and p^ax are the maximal velocity and the maximal density of pedestrians respectively. Throughout this paper 
the LWR flux will be used. Next p: R+ ^ R+ prescribes the maximal flow allowed through a bottleneck located at x = 0 
as a function of the weighted average density in a left neighbourhood of the bottleneck and w: ^ K+ is the weight 

function used to average the density. 

Finally in (fT^ . p(f, 0-) denotes the left measure theoretic trace along the constraint, implicitly defined by 


Y ^+oo ^0 

lim- I I |p(f, x)-p(f, 0-)| 0(f, x) dxdf = 0 

efo e Jo J-£ 


for all e C“(M^ M). 
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The right measure theoretic trace, p(t, 0+), is defined analogously. 


at locations with reduced 
I3> 11 1 present 


In the last few decades, the study of the pedestrian behaviour through bottlenecks, namel 
capacity, such as doors, stairs or narrowings, drawn a considerable attention. The papers ll5| 
results of empirical experiments. However, for safety reasons, experiments reproducing extremal conditions such as 
evacuation and stampede are not available. In fact, the unique experimental study of a crowd disaster is proposed in m. 
The available data show that the capacity of the bottleneck (i.e. the maximum number of pedestrians that can flow through 
the bottleneck in a given time interval) can drop when high-density conditions occur upstream of the bottleneck. This 
phenomenon is called capacity drop and can lead to extremely serious consequences in escape situations. In fact, the 
crowd pressure before an exit can reach very high values, the efficiency of the exit dramatically reduces and accidents 
become more probable due to the overcrowding and the increase of the evacuation time (i.e. the temporal gap between the 
times in which the first and the last pedestrian pass through the bottleneck). A linked phenomenon is the so-called Faster 
Is Slower (FIS) effect, first described in iQ. FIS effect refers to the jamming and clogging at the bottlenecks, that result 
in an increase of the evacuation time when the degree of hurry of a crowd is high. We recall that the capacity drop and the 
FIS effect are both experimentally reproduced in la 1 141] . A further related (partly counter-intuitive) phenomenon is the 


so-called Braess ’ paradox for pedestrian flows ||15|] . It is well known that placing a small obstacle before an exit door can 
mitigate the inter-pedestrian pressure and, under particular circumstances, it reduces the evacuation time by improving 
the outflow of people. 

Note that as it happens for any first order model, see for instance ifl^ Part III] and the references therein, ADR can not 
explain the capacity drop and collective behaviours at the bottlenecks. Therefore one of the difficulties we have to face is 
that the constraint p has to be deduced together with the fundamental diagram from the empirical observations. 

The aim of this paper is to validate ADR by performing simulations in order to show the ability of the model to 
reproduce the main effects described above and related to capacity drop that are FIS and Braess’ paradox. To this end 
we propose a numerical scheme for the model and prove its convergence. The scheme is obtained by adapting the local 
constrained finite volume method introduced in 10 to the non-local case considered in ADR, using a splitting strategy. 

The paper is organized as follows. In Section we briefly recall the main theoretical results for ADR. In Section [3] 
we introduce the numerical scheme, prove its convergence and validate it with an explicit solution obtained in 111]. In 
Section |4] we perform simulations to show that ADR is able to reproduce the Braess’ paradox and the FIS effect. In 
Subsection l4.3l we combine local and non-local constraints to model a slow zone placed before the exit. Conclusions and 
perspectives are outlined in Section|5] 


2. Well-posedness for the ADR model 

Existence, uniqueness and stability for the general Cauchy problem ([T]) are established in |lt] under the following 
assumptions: 

(F) / belongs to Lip ([0,R]; [0, -HooQ and is supposed to be bell-shaped, that is /(O) = 0 = f(R) and there exists 
cr e ]0,R[ such that f'(p) (cr -p) > 0 for a.e. p e [0,R]. 

(W) w belongs to L“’(]R_;M+), is an increasing map, ||w||l 1 (k j = 1 and there exists iw > 0 such that w{x) - 0 for 
any .r < -i,^- 

(P) p belongs to Lip ([0, R]; ]0, /(cr)]) and is a non-increasing map. 

The regularity w 6 L°°(K_; R+) is the minimal requirement needed in order to prove existence and uniqueness of ([TJ. 
In this paper, we shall consider continuous w. 


The existence of solutions for the Riemann problem for ([TJ is proved in ll8|] for piecewise constant p. However, such 
hypothesis on p is not sufficient to ensure uniqueness of solutions, unless the flux / and the efficiency p satisfy a simple 
geometric condition, see MlSll for details. In the present paper, we consider either continuous nonlinear p or a piecewise 
constant p that satisfies such geometric condition. 

The definition of entropy solution for a Cauchy problem (ITaJ . (ITbJ with a fixed a priori time dependent constraint 
condition 


/(p(f,0+)) < q{t) 


t 6 . 


( 2 ) 


was introduced in E Definition 3.2] and then reformulated in E Definition 2.1], see also E Proposition 2.6] and | 
Definition 2.2]. Such definitions are obtained by adding a term that accounts for the constraint in the classical definition 
of entropy solution given by Kruzkov in |21, Definition 1]. The definition of entropy solution given in |ll Definition 2.1] 
is obtained by extending these definitions to the framework of non-local constraints. 

The following theorem on existence, uniqueness and stability of entropy solutions of the constrained Cauchy prob¬ 
lem ([TJ is achieved under the hypotheses (F), (W) and (P). 
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Theorem 2.1 (Theorem 3.1 in |[ll]). Let (F), (W), (P) hold. Then, for any initial datum p e L“(R; [0,/?]), the Cauchy 
problem ([TJ admits a unique entropy solution p. Moreover, if p' — p'(t,x) is the entropy solution corresponding to the 
initial datum p' e L“(R; [0, R]), then for all T > 0 and L > i^, the following inequality holds 

where M = Lip(/) and C = 2Lip(/?)||w||L«.(R ). 

The total variation of the solution may in general increase due to the presence of the constraint. In (ll] the authors 
provide an invariant domain D c L,^ (R; [0,/?]) such that if p belongs to D, then one obtains a Lipschitz estimate with 
respect to time of the norm and an a priori estimate of the total variation of 


'T(p) = sign(p - o-)[/(cr) - /(p)] 



dr. 


3. Numerical method for approximation of ADR 

In this section we describe the numerical scheme based on finite volume method that we use to solve ([1]). Then we 
prove the convergence of our scheme and validate it by comparison with an explicit solution of ([1]). In what follows, we 
assume that (F), (W) and (P) hold. 


3.1. Non-local constrained finite volume method 

Let Ajc and Af be the constant space and time steps respectively. We define the points Xy+ 1/2 = y'Ajc, the cells Kj - 
{xj-ii 2 ,Xj+ii 2 { and the cell centers Xj - (j - l/2)Ax for j e Z. We define the time discretization f” = nAt. We introduce 
the index jc such that Xj^+ 1/2 is the location of the constraint (a door or an obstacle). For n e N and j e Z, we denote by 
p" the approximation of the average of p{f \ ■) on the cell Kj, namely 


o_^ 


p(x) Ax 


and 


\ r^^j+ip 


p(f, x) Ax 


if n > 0. 


We recall that for the classical conservation law (fTali-(fTb]). a standard finite volume method can be written into the form 

(4) 

where ~ ^ {p'j’P%i) ^ nionotone, consistent numerical flux, that is, F satisfies the following assumptions: 

• F is Lipschitz continuous from [0, Rf to R with Lipschitz constant Lip(F), 

• F(a, a) — f(a) for any a e [0, R], 

• (a, b) e [0, R^ ^ F(a, b) eM.is non-decreasing with respect to a and non-increasing with respect to b. 

We also recall that in iB the numerical flux for the time dependent constraint (|2]i is modified as follow in order to take 
into account the constraint condition 


<x-« 

^ 1 + 1/2 


min{F(p],p”_^i),^''} if j^jc. 


(5) 


where q'‘ is an approximation of q(f ). In the present paper, when dealing with a Cauchy problem subject to a non-local 
constraint of the form (fTcl) we will use the approximation 


q =p 


Ax ) w(xj)p" 


j^jc 


( 6 ) 


Roughly speaking 

• we apply the numerical scheme Q for the problem (fTal)- (fTb]) . 

• we apply the numerical scheme (m-Q for the problem (fTal)- (llbl i-(l2Ti. 

• we apply the numerical scheme (|4|-(|5]l-® for the problem ([1]). 
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3.2. Convergence of the scheme 

Let us introduce the finite volume approximate solution pA defined by 

PA(f, x) - p" for X 6 Kj and t e [f, [, (7) 

where the sequence is obtained by the numerical scheme (Si-©. Analogously, we also define the approximate 

constraint function 


qk(.t) = q’' 


for t e [f, [. 


( 8 ) 


First, we prove a discrete stability estimate valid for any domain 2 = [0, L] x K with T > 0, for the scheme ©-© 
applied to problem (fTali- dlbl i-dZt. This estimate can be seen as the equivalent, in this framework, of the stability result 
established in O Proposition 2.10]. 

Proposition 3.1. Letp be in L°°(]R; [0,/?]) and q^, q^ be piecewise constant functions of the form ©. If Pa and pA are 
the approximate solutions of ([Ial)-([Tb]i-© corresponding, respectively, to q^ and q^ and constructed by applying the 
scheme @-©, then we have 

IIPa “PaIIl*(0 ^ 27’II^a - §AllLi([0,r])- 

Proof. For notational simplicity, let N -\T /AfJ. Let us also introduce (pp^ez.neN defined by, 

~pT' ^P"j-Tx - ^J-m) - for any 7 e Z, n 6 N, 


where is defined by 


qrn 

^7+1/2 


min {f (p^pPi) . §"} if j = jc- 


Then using the definitions of (p'fjez.neN and (p''-)jez,neN, we have for any n = 1,..., A, 


p^^p" if 

j i [jc, 7c + 1) 

and 



+ min(F(p”;',c’-',).r'|). 

pi,i-Pl.t = ^ 


which implies that 


\pl-pl\^^\r'-r'\. 

K., -P1..I ^ s W 

Therefore we deduce that, for any n - 1 ,..., A, 



n—l I 




I 


(9) 


Besides, observe that the modification of the numerical flux at the interface Xj^+ 1/2 introduced in © does not affect the 
monotonicity of the scheme ©-© (see iB Proposition 4.2]). Therefore, for any n - 1,..., A, we have 




(10) 


Hence thanks to © and (fTOl) . we can write 


z K- -pj\ ^ z K -pj ^+z -^‘1 - z ^p'j-p'j ^= 


At 


Ax 


Ax 


f - fl 


Then an induction argument shows that for any n - 1,..., A, 




jeZ 




II^A - §AllL>([0,f"])- 


4 


In conclusion, we find that 


IIPa - Pa\\l'(Q) - ^ ^ Ip" ~ p’jl ^ 2||^a - §AllLi([0,ri) ^ ^ 27’||q'A - §AllL>([0,r]) 

n=l yeZ n=l 

and this ends the proof. 

Let us now notice that as in ini Proposition 4.2], under the CFL condition 

At 1 
Ltp(F) - , -, 


□ 


for a.e. (f, x) e Q. 


( 11 ) 

( 12 ) 


we have the L“ stability of the scheme (|4]i-(|5]l-(l6]l that is 

0 < p^{t, x) < R 
This stability result allows to prove the statement below. 

Proposition 3.2. Let be defined by (|6])-([8]l. Then under the CFL condition dill) , for any T > 0, there exists C > 0 only 
depending on T, f, F, p, w and R such that: 


l?AlBV([0,r]) - 

Proof. Let N - \ T /AfJ and be an integer such that supp(vi.') c U Kj. Then for any n = 0,..., - 1, we have 

h<i<ic 


(13) 




Ax ^ w{xj)pf^ -p Ax w(xj)p"j 

h<j<jc / jw<j<jc 


< AxLip(p) 

j.,<j<jc 

Now, using a summation by part, we have 




At Lip(p) 




h<j<ic 


Z ('^>1/2 “ ^"-1/2) = ^(■*l)'^"+i /2 - - Z (>^(■^2+1) - ^(-^f)) 


2+1/2- 


h<j<jc 

Then, it follows that 


L<j<jc-i 


1^"+' - ^"1 < Af Lip(p) ||w||l~(r_;R) Z I^Li/2l- 

L-l<j<jc 

Now, from Q, for any j eh we have the estimate 

^ k(PpP}ri)| ^ ^(PpPiri) - ^(PpP’Pl + |/(Ppl ^ Lip(F) |p;Vi - p"\ + Lip(/) Ip'JI < R (Lip(F) + Lip(/)) . 
Hence we deduce that 


N-l 


I?aIbv([0,7']) - Z 1^"^' ~ ^"1 ~ 

rt =0 


where C = (jc - i„+2)TR Lip(p) | 


(Lip(F) + Lip(/)). 


□ 


We are now in a position to prove a convergence result for the scheme dH)-©-®. 

Theorem 3.1. Under the CFL condition dill) , the constrainted finite volume scheme ©-©-© converges in L^{Q) to the 
unique entropy solution to ©. 

Proof Let (pA,qA) be constructed by the scheme ©-©-©. Proposition 13.21 and Helly’s lemma give the existence of a 
subsequence, still denoted qA and a constraint function q e L°°([0, T]) such that qA converges to q strongly in L^([0, T]) 
as At —» 0. Let p 6 L”(]R+ x R; [0, R]) be the unique entropy solution to (fTal)- dlbb -(l2l) associated to q. It remains to prove 
that the subsequence pA converges to p strongly in L^(2) as At, Ax —> 0. The uniqueness of the entropy solution to (fTab - 
(ITbl) -(l2l) will then imply that the full sequence pA converges to p and, as a consequence, the full sequence qA converges to 
P(L- w(x)p(f,x)dx). 

Let qA be a piecewise constant approximation of q such that qA converges to q strongly in L^([0, T]). Furthermore, we 
also introduce pA constructed by the scheme ©-© and associated to qA- Now we have 


lip “ paIIli( 0 5 lip “ paIIli(g) + IIpa “ paIIlho)- 

But, thanks to IB Theorem 4.9], under the CFL condition (fTTT) . ||p -PaIIl^g) tends to 0 as At, Ax —> 0. Furthermore, 
thanks to Proposition l3.ll we have 

IIPa - PaIIlI(G) -2 7’ Ilg'A - §aIIl1([0,7']) 

which also shows that ||pA - PaIIlhg) l^tids to 0 as At, Ax —» 0. □ 
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5.5. Validation of the numerical scheme 



Figure 1: The functions [p h-> /(p)] and [f i-> p(f)] as in Section [331 


We propose here to validate the numerical scheme (|4]i-(|5])-(|6]l using the Godounov numerical flux (see e.g. 1122112311 1 
which will be used in the remaining of this paper: 


F{a, b) 


( min/ 

1 la.b] 

1 maxf 

[ M 


if a < b, 
if a > b. 


We consider the explicit solution to ([TJ constructed in (ll Section 6] by applying the wave front tracking algorithm. The 
set up for the simulation is as follows. Consider the domain of computation [-6,1], take a normalized flux f(p) = p(l -p) 
(namely the maximal velocity and the maximal density are assumed to be equal to one) and a linear weight function 
w{x) = 2(1 + x)x[^i,o]{x). Assume a uniform distribution of maximal density in [xa, xb\ at time f = 0, namely p - X[xa.xbV 
The efficiency of the exit, p, see Figure [T] is of the form 


PiO = 


The explicit solution p corresponding to the values 

po = 0.21, Pi =0.168, p2= 0.021, 


PO 

if0</</i. 



Pi 

if/l </</2, 



P2 

if/2</< 1. 




~ 0.566, Xa = -5.75, 

xb - -2, 

ft -0.731, 


is represented in Figure |2l The above choices for the flux / and the efficiency p ensure that the solution to each Riemann 
problem is unique, see ifisll . We defer to IH Section 6] for the details of the construction of the solution p and its physical 
interpretation. 

A qualitative comparison between the numerically computed solution x i-> pa(6 x) and the explicitly computed solution 
X p(f, x) at different fixed times t is in Figure [3] We observe good agreements between x i-> p(f, x) and x pa(6 x). 
The parameters for the numerically computed solution are Ax = 3.5 x 10 and At - 7 x 10“^. 

A convergence analysis is also performed for this test. We introduce the relative L^-error for the density p, at a given time 
f, defined by 


E" - 


Y^\p(f\xj)-p’}\ I ^|p(f",Xj)| 


In Table [U we computed the relative -errors for different numbers of space cells at the fixed time f = 10. We deduce 
that the order of convergence is approximatively 0.906. As in d, we observe that the modification (|5]l of the numerical 
flux does not affect the accuracy of the scheme. 


4. Numerical simulations 

This section is devoted to the phenomenological description of some collective effects in crowd dynamics related to 
capacity drop, namely the Braess’ paradox and the Faster Is Slower (FIS) effect. 


4.1. Faster is Slower effect 

The FIS effect was first described in |[T3 , 24 ] in the context of the room evacuation problem. The authors studied the 
evolution of the evacuation time as a function of the maximal velocity reached by the pedestrians, and they shown that 
there exists an optimal velocity for which the evacuation time attains a minimum. Therefore, any acceleration beyond the 
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(a) The solution in the (/, jc,p)-coordinates. 


(b) The solution in the (x, 0-coordinates. 




- 0.2 


(c) The solution in the (?, x,p)-coordinates for 0 < t< 15. 


0.0 

(d) The solution in the (?, .j:, p)-coordinates for 85 < t < 87.5. 


- 0.4 


Figure 2: Representation of the solution constructed in Q], Section 6] and described in Subsection l3.3l 
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(a) pa{0,x) 


(b) pa(1,-c) 


(c) pa(7.325,x) 



(d)pA(10,;f) 


l.On -1- 

Q g- 

Q g- 

- 

Q 2- 

O.oU_^^^_ I_^^_ 

-6 -5 -4 -3 - 2-1 0 1 

(e) p(0,;f) 



(f) P(hx) 



(g) p(1.325,x) 



(h) p(10,x) 



(i) pa(1 1.939, X) 



(j) Pa(85.2,x) 



(k) pa(85.5526,x) 


(1) pa(87.4981,x) 



-6 -5 -4 -3 -2 -1 

(m) p(l 1.939, x) 


-6 -5 -4 -3 - 2-10 1 

(n) p(85.2,x) 


-5 -4 - 3 - 2-1 0 1 

(o) p(85.5526,x) 


- 6-5 -4 -3 - 2-1 0 1 

(p) p(87.4981,x) 


Figure 3: With reference to Subsection l3.3l The numerically computed solution x i-> jc) and the explicitly computed solution x p{t, x) at different 
fixed times t. 
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Number of cells 

L^-error 

625 

9.6843 X 10-^ 

1250 

6.2514 X 10-^ 

2500 

3.4143 X 10-^ 

5000 

1.3172 X 10-^ 

10000 

1.03 X 10-^ 

20000 

4.2544 X 10-"^ 

Order 

0.906 


Table 1: Relative -error at time t = 10. 


optimal velocity worses the evacuation time. Following the studies above, the curve representing the evacuation time as a 
function of the average velocity takes a characteristic shape 1124 Figure Ic]. 

The first numerical tests we performed aim to verify if such shape is obtained starting from the ADR model. To this end, 
we consider the corridor modeled by the segment [-6,1], with an exit at x = 0. We consider the flux f{p) - p Vmax (1 - p) 
where Vmax is the maximal velocity of the pedestrians and the maximal density is equal to one. We use the same weight 
function as for the validation of the scheme, w{x) - 2(1 -i- x)x[-ifl](x) and, the same initial density, p - [_ 5 . 75 __ 2 ]■ The 

efficiency of the exit p is now given by the following continuous function 


P(0 


Po 

iPQ - Pl)^ + Pl^l - P0^2 


^1 - ft 


Pi 


if0<f <ft, 
if ft <^<ft, 
ifft<^<l, 


(14) 


where 


Po - 0.24, Pi - 0.05, ft = 0.5, ft = 0.9. 


The space and time steps are fixed to Ax = 5 x 10 ^ and Af = 5 x 10 In Figure|4]are plotted the flux / corresponding 
to the maximal velocity Vmax = 1 and the above efficiency of the exit. 



P. f 


Figure 4: The normalized flux p —» f(p) and the constraint ^ —» p(^) defined in <141 . 

Figure [^represents the evacuation time as a function of the maximal velocity as Vmax varies in the interval [0.1,5]. 
As we can observe, the general shape described above is recovered. The numerical minimal evacuation time is 19.007 
and is obtained for Vmax = 1 ■ 

In addition, we reported in Figure [6] the density at the exit as a function of time for different values of the maximal 
velocity Vmax around the optimal one. We notice that the maximal density at the exit and the time length where the density 
is maximal increase with the velocity. This expresses the jamming at the exit that leads to the FIS effect. 

Then we performed some series of tests to see how the general shape obtained in Figure [5] changes with respect to 
variations of the parameters of the model. In Figure |8| (a), we show this variation when we consider different initial 
densities, namely, p, pi and ft with pi(x) = 0.8;^f[_5.75__2] and ft(x) = 0.6;if[_5.75__2]- The general shape of the curves is 
conserved. We observe that the evacuation time increases with the initial amount of pedestrians while the optimal velocity 
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Figure 5: With reference to Subsection l4. II Evacuation time as a function of the velocity v„,ax- 




t t 

(a) p i-> Pa(0, 0 for velocities <1. (b) p i-> Pa(0, t) for velocities v^ax S 1. 


Figure 6: With reference to Subsec tion l4. II Densities at the exit as a function of time for different velocities. 


decreases as the initial amount of pedestrians increases. The minimal evacuation time and the corresponding optimal 
maximal velocity are 12.259 and 1.07 forp 2 and 15.691 and 1.03 forpi. 

Next we explore the case where the efficiency of the exit varies. We consider the function p defined in (fT4l) and the 
modification pp such that pp{^) - p(J30- In Figure [T] we plotted the functions p, pp for y6 = 0.8 and p - 0.9. Then, 
in Figure [8] (b) are plotted the evacuation time curves corresponding to these three efficiencies of the exit. As minimum 
evacuation times, we obtain 18.586 and 18.827 for^S = 0.8, 0.9 respectively. As expected, the minimal evacuation time 
increases with lower efficiency of the exit. The corresponding velocities are approximatively 1.06 and 1.02 respectively. 

Finally, we change the location of the initial density. In addition to the corridor [-6,1], we consider two other corridors 
modeled by the segments [-12,1] and [-20,1]. In these two corridors we take as initial densities pz{x) - ;t'[-ii. 75 -g] 
Pa{x) - ;t'[-i 9 . 75 - 16 ] respectively. We have reported the obtained evacuation time curves in Figure[8](c). As expected, the 
minimal evacuation time increases with the distance between the exit and the initial density location. 

4.2. Braess ’ paradox 

The presence of obstacles, such as columns upstream from the exit, may prevent the crowd density from reaching 
dangerous values and may actually help to minimize the evacuation time, since in a moderate density regime the full 
capacity of the exit can be exploited. From a microscopic point of view, the decrease of the evacuation time may seem 
unexpected, as some of the pedestrians are forced to chose a longer path to reach the exit. 

The ADR model is able to reproduce the Braess’ paradox for pedestrians, as we show in the following simulations. We 
consider, as in the previous subsection, the corridor modeled by the segment [-6,1] with an exit at x = 0. We compute 
the solution corresponding to the flux f(p) = p(l - p), the initial density p(x) = efficiency of the exit p 

of the form (I14l i with the parameters 


Po - 0.21, 


Fi =0.1, 


10 


= 0.566, 


6 = 0.731 


















































0.25 


0.2 


0.15 


0.1 


0.05 



Figure 7: With reference to Subsection l4.ll The efficiencies ^ for yS = 0.8,0.9,1. 


and the same weight function w{x) - 2(1 + x)x[-i,o](x). The space and time steps are hxed to Ajc = 5 x 10“^ and 
Af = 5 X 10 Without any obstacle, the numerical evacuation time is 29.496. In these following simulations we place an 
obstacle atx = d, with -2 < d < 0. The obstacle reduces the capacity of the corridor and can be seen as a door, which we 
assume larger than the one at x = 0. Following these ideas we dehne an efficiency function pdi^) = 1.15 p{^) and a weight 
function Wd(x) - 2(x - d + y)X[d-\,d]{x) associated to the obstacle. 

In Figurel^we have reported the evolution of the evacuation time when the position of the obstacle varies in the interval 
[-1.9, -0.01] with a step of 0.01. We observe that for -1.8 < d < -1.72, the evacuation time is lower than in the absence 
of the obstacle. The optimal position of the obstacle is obtained for d - -1.72 and the corresponding evacuation time 
is 24.246. We compare in Figure [10] hve snapshots of the solution without obstacle and the solutions with an obstacle 
placed ai d - -1.72 and d - -1.85. This latter location corresponds to a case where the evacuation time is greater than 
the one without an obstacle. In these snapshots, we see that the obstacle placed al d - -1.85 becomes congested very 
soon. This is due to the fact that the obstacle is too close to the location of the initial density. When the obstacle is placed 
aid - -1.72, it delays the congestion at the exit. 

4.3. Zone of low velocity 

In this section, we perform a series of simulations where the obstacle introduced in Subsection 14.21 is now replaced by 
a zone where the velocity of pedestrians is lower than elsewhere in the domain. The effect we want to observe here is 
similar to the one we see in Braess’ Paradox. Namely we prevent an high concentration of pedestrians in front of the 
exit by constraining their flow in an upstream portion of the corridor. In this case however the constraint is local, as the 
maximal value allowed for the flow only depends on the position in the corridor. 

We consider again the corridor modeled by the segment [-6,1] with an exit at x = 0. The efficiency of the exit and the 
initial density are the same as in the previous subsection. Assume that the slow zone is of size one and is centred atx-d, 
where -1.9 < < 0. Define the following function 


k(x) 


1 

-2(x - d) 
2(x - d) 

1 


if X < c/ - 0.5, 
if d -0.5 < X < d, 

if d < X < d + 0.5, 
if x>d + 0.5, 


(15) 


and the following velocity v(x,p) = [T + (1 - T) ^(x)] Vmax (1 - p), where 4 e [0,1] and Vmax > 1 is the maximal velocity. 
With such velocity, the maximal velocity of pedestrians decreases in the interval {d - 0.5, d}, reaching its minimal value 
4 Vmax atx-d. Then the velocity increases in the interval {d, d + 0.5] reaching the maximum value Vmax, that corresponds 
to the maximal velocity away from the slow zone. Finally we consider the flux /(x,p) = p v(x,p) and the space and time 
steps are fixed to Ax = 5 x 10“^ and At - 5 x 10“"^. 

Figure [TT] (a) shows the evolution of the evacuation time as a function of the parameter A varying in the interval [0.1,1] 
when the center of the slow zone is fixed at d - -1.5. We observe that the optimal minimal velocity in the slow zone is 
for A - 0.88 and the corresponding evacuation time is 20.945. Recalling that without the slow zone the evacuation time is 
29.496, we see that the introduction of the slow zone allows to reduce the evacuation time. In Figure [TT](b), we show the 
evolution of the evacuation time when varying the center of the slow zone d in the interval [-1.9,0] and when the minimal 
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(a) Evacuation time as a function of v^ax for different amounts of initial 
densities. 



(b) Evacuation time as a function of v^ax for different efficiencies of the 
exit. 



Vmax 

(c) Evacuation time as a function of Vmax for different locations of the 
initial density. 


Eigure 8: With reference to Subsection 14. II Evacuation time as a function of v^ax for different parameters of the model. 


and the maximal velocities are fixed and correspond to d = 0.88 and Vmax = 1- We observe here that, unlike in the Braess 
paradox tests case, the evacuation time does not depend on the location of the slow zone, except when this latter is close 
enough to the exit. Indeed, when the slow zone gets too close to the exit, the evacuation time grows. This is due to the 
fact that pedestrians do not have time to speed up before reaching the exit. 

Fix now d — -1.5 and A - 0.88 and assume that Vmax varies in the interval [0.1,5]. The evolution of the evacuation 
time as a function of Vmax is reported in Figure[TT](c). We observe that we get the characteristic shape already obtained in 
the FIS effect. 

Finally we present in Figure[T2]five snapshots for three different solutions. The first two solutions are the ones computed 
in Subsection l4.2l without obstacle and with an obstacle located atd = -1.72 respectively. The third solution is computed 
with a zone of low velocity centered atd = -1.72, d = 0.88 and Vmax = 1. In order to have a good resolution of this third 
solution, the space and time steps where fixed to Ax = 3.5 x 10“^ and Af = 7 x 10“^.We note that in the case where a zone 
of low velocity is placed in the domain, we do not see the capacity drop, as the density of pedestrians never attains very 
high values in the region next to the exit. 


5. Conclusions 

Qualitative features that are characteristic of pedestrians’ macroscopic behaviour at bottlenecks (Faster is Slower, 
Braess’ paradox) are reproduced in the setting of the simple scalar model with non-local point constraint introduced in (ll] . 
These effects are shown to be persistent for large intervals of values of parameters. The validation is done by means of a 
simple and robust time-explicit splitting finite volume scheme which is proved to be convergent, with experimental rate 
close to one. 

The results presented in this paper allow to consider more complex models. Indeed, as ADR is a first order model, 
it is not able to capture more complicated effects related to crowd dynamics. Typically, ADR fails to reproduce the 
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Figure 9: With reference to Subsection l4.2l Evacuation time as a function of the position of the obstacle. 


amplification of small perturbations. This leads to consider second order model such as the model proposed by Aw, 
Rascle and Zhang ll25[ 12611 in the framework of vehicular traffic. 

Another extension of this work is to consider the ADR model with constraints that are non-local in time. Such con¬ 


straints allow to tackle optimal management problems in the spirit of Il27i 12811 . 

Finally, this work can also be extended to two-dimensional models where experimental validations may be possible. 
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Figure 10: With reference to Subsection l4.2l Braess paradox simulations: density profiles at times t = I (first line), t = 1 (second line), f = 15 (third 
line), f = 19 (fourth line) and t = 24.246 (last line). 
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(a) Evacuation time as a function of A. 
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(b) Evacuation time as a function of d. 
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(c) Evacuation time as a function of Vmax- 


Figure 11: With reference to Subection l4.3l Evacuation time as a function of different parameters of the model. 
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Zone of low velocity centered at 
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Figure 12: With reference to Subsection 14.31 Braess’ paradox and zone of low velocity simulations: density profiles at times f = 1 (first line), t = 1 
(second line), f = 15 (third line), f = 19 (fourth line) and t = 20.945 (last line). 
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